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Abstract 

We study a solution of long polyanions (PA) with shorter polycations (PC) and focus on the role 
of Coulomb interaction. A good example is solutions of DNA and PC which are widely studied 
for gene therapy. In the solution, each PA attracts many PCs to form a complex. When the ratio 
of total charges of PA and PC in the solution, x, equals to 1, complexes are neutral and they 
condense in a macroscopic drop. When x is far away from 1, complexes are strongly charged. The 
Coulomb repulsion is large and free complexes are stable. As x approaches to 1, PCs attached to PA 
disproportionate themselves in two competing ways. One way is inter-complex disproportionation, 
in which PCs make some complexes neutral and therefore condensed in a macroscopic drop while 
other complexes become even stronger charged and stay free. The other way is intra-complex 
disproportionation, in which PCs make one end of a complex neutral and condensed in a small 
droplet while the rest of the complex forms a strongly charged tail. Thus each complex becomes 
a "tadpole". These two ways can also combine together to give even lower free energy. We get a 
phase diagram of PA-PC solution in a plane of x and inverse screening radius of the monovalent 
salt, which includes phases or phase coexistence with both kinds of disproportionation. 

PACS numbers: 61.25.Hq, 82.35.Rs, 87.14.Gg, 87.15.Nn 
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I. INTRODUCTION 



Condensation in solution of two oppositely charged polyelectrolytes (PE) is an important 
phenomenon in biology and chemical engineering. One of the most interesting applications 
is DNA condensation by polycations (PC), which is widely^used in gene therapy research. 
A good example is condensation of DNA with poly-lysine 1]. Complexation of DNA with 
PCs can invert the charge of bare DNA and help DNA to penetrate negatively charged cell 
membrane. At the same time, adsorbed PC in complexes or their condensate may protect 
DNA from digestion by enzymes inside the cell Tremendous amount of experimental 
works have been done in this area jll Q, 0, S, aT" 



(a) 



(b) 







FIG. 1: Objects appearing in a solution of PA and PC (left column) and their symbols used in 
Figs. [21 HI 1^ and IIUI (right column). The long polymer is PA and the short polymer is PC. (a) a 
single PC. (b) negative PA-PCs complex, (c) positive PA-PCs complex, (d) condensate of almost 
neutral complexes, (e) tadpole made of one PA-PCs complex. Here only the case of positive tail 
is shown. The tail can be negative, too. 

In this paper, motivated by DNA-PC condensation, we study the equilibrium state of a 
solution of polyanions (PA) and polycations (PC) in the presence of monovalent salt and 
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focus on the role of Coulomb interaction. We assume that both PA and PC are so long that 
at room temperature T their translational entropy can be ignored in comparison with their 
Coulomb energy. We are particularly interested in the case when PA is much longer than PC 
and many PCs are needed to neutralize one PA. In Fig. Qwe list objects which can appear in 
such a solution. Each PA attracts many PCs to form a PA-PCs complex (Fig.^,c). Neutral 
complexes can further condense in a liquid drop (Fig. ^). One complex can form a neutral 
head and a charged tail to become a tadpole (Fig. [T^). And it is possible to have excessive 
free PCs (Fig. ^). When and where these objects exist or co-exist with each other depends 
on two dimensionless parameters: the ratio of total charges of PC and PA in the solution, 
X, and b/vg, where b is the size of the monomer and is the Debye-Hiickel screening radius 
provided by monovalent salt. The main result of this paper is the phase diagram in a plane 
of X and b/vs as shown in Fig. |21 We discover a new phase of "tadpoles" originating from 
the polymer nature of the objects. We also present a first theory of broadening of the phase 
of a single drop with decreasing (curves x^lvg) and 2:4 (r^) in Fig. |2)). 

Up to now, there was no complete theory of phase diagram for such systems. Previously, 
the complexation of oppositely charged polyelectrolytes was studied in a symmetric system 
in which the length, concentration and linear charge density of PA and PC are the same 
It was shown that even in the absence of monovalent salt, strongly charged PA and PC form 
a single macroscopic drop of neutral dense liquid, which separates from water. It corresponds 
to the phase at x = 1 in our phase diagram (Fig. [2|). On the other harid, the phase diagram 
of a solution of DNA and short polyamines was studied in Refs. in which 

translational entropy of polyamines plays a very important role. We postpone discussion 
of the role of translational entropy till the end of this paper (Sec. IVI|) . In our previous 
works 3; M\, phase diagrams have been discussed for other systems. In Ref. [isl, solution 
of very long PA with positive spheres was considered. This system is similar to chromatin 
in the sense that each PA binds many spheres making a long necklace. We also discussed 
the phase diagram of a system of oppositely charged spheres in strongly asymmetric case 
when each say negative sphere complexes with many positive ones jl4|. Many features of 
the phase diagram in Fig. |21 are also applicable to these systems and we will return to them 
in the conclusion. 

A recent paper of us studied the exactly same topic 15|- But it missed two phase 
coexistence regions at ^ 00 (comparing Fig. 2 of Ref. with Fig. |21in this paper). The 
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FIG. 2: Typical phase diagram of a solution of PA and PC. The horizontal axis x is the ratio of 
total charges of PC and PA in the solution. The vertical axis b/rg is the ratio of the length of a 
monomer size of the PA molecule to the Debye-Hiickel screening radius r^. Symbols are explained 
in Fig. n 

main purpose of the present paper is to make corrections and to discuss the complete phase 
diagram we have right now. 

It is enough to look at x < 1 side of Fig. |21to see the main feature of the phase diagram. 
It contains three phases: the phase of a single drop for x > x^, the phase of free complexes 
for X < xi, and the "tadpole" phase for X2 < x < X3. Between these phases, there are 
three regions of phase coexistence. This phase diagram is very much similar to the phase 
diagram of water in temperature and volume coordinates jl^, when r^/fo is considered as 
temperature and x as volume. The three phases mentioned above are like gas, solid and 
liquid respectively. Essentially, this analogy originates from the Gibbs' phase rule , which 
is crucial to determine our phase diagram (see Sec. IIII|) . 

Let us now try to understand the physics of this phase diagram. We start from the 
horizontal axis (r^ 00) and first focus on x < 1 side. In the solution, each PA adsorbs 
many PCs to form a complex. When x <^ 1, the number of PC is not enough to neutralize 
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all PAs and each PA-PCs complex is strongly negatively charged (Fig. QJd). The Coulomb 
repulsion between complexes is huge and all complexes stay free, or in other words, colloidal 
solution of complexes is stable (see ranges < x < Xi in Fig. |2I). 

When a; = 1, each PA-PCs complex is neutral and there is no Coulomb repulsion between 
them. They all condense to form a macroscopic strong correlated liquid drop (see Fig. 
Due to the orderly arrangement of positive and negative charges in the drop, a certain 
amount of short-range correlation energy is gained (Fig. E)). We define e < as the energy 
gain of a neutral complex in the macroscopic drop. Note that monomers on the surface 
of the drop can not gain as much energy as monomers inside. This defines the surface 
energy of the drop which plays an important role in the competition between two kinds of 
disproportionations (see the next paragraph). 




FIG. 3: A schematic illustration of short-range attraction between neutral complexes in a condensed 
liquid drop using the model of strongly charged PA and PC. A portion of two complexes in the 
liquid drop is shown. PA and PC charges are shown by black and white dots correspondingly. 
The dashed lines show two complexes sitting in parallel planes. The complexes attract each other 
because charges of the same sign are farther away than charges of the opposite sign. 

In vicinity of a; = 1, the long-range Coulomb repulsion between charged complexes com- 
petes with the short-range attraction due to correlations. As x increases, condensation starts 
at X = Xi. To minimize the free energy, PCs can be redistributed among complexes so that 
a portion of complexes are neutral and therefore condensed in a macroscopic drop, while 
the rest of complexes become stronger charged and stay free. This is called inter- complex 

n rn 

disproportionation \2\ or partial condensation [1^ (see Fig. H^). It is essentially a coexis- 
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tence of the two phases: the single drop phase and the free complexes phase. On the other 
hand, PCs can also disproportionate themselves within each complex, which we call intra- 
complex disproportionation (Fig. EJd). PCs move closer to one end of PA molecule, making 
one part of a PA molecule neutral and condensed in a droplet, while the other part is even 
stronger charged and not condensed. Unlike inter-complex disproportionation, this gives a 
new "tadpole" phase, as far as L is finite, where L is the length of a free PA-PCs complex 
(see Fig. 121). It is possible to combine the two ways of disproportionation to accomplish even 
lower free energies (see Fig. Eh,d). As a result, at — * oo, we have the sequence of these 
phases as shown on the horizontal axis of Fig. |21 

What happens at x > 1 side is almost the same, except that here the number of PC is 
larger than necessary to neutralize all PAs and each complex is strongly positively charged 



(charge inverted) ( 
sion at a; > 1 side 



^). Let us briefly remind the nontrivial mechanism of charge inver- 
17[ . We illustrate it in Fig. El for the model of strongly charged flexible 



PA and PC, in which the distance between charges, b, is the same for PA and PC molecules 
and is of the order of Bjerrum length Ib = 7A [e'^/Dls = ksT, Z) = 80 is the dielectric 
constant of water). When a new PC molecule arrives at a neutral PA-PCs complex, all PCs 
in the complex can rearrange themselves so that the charge of this excessive PC is smeared 
in the whole complex and the Coulomb self-energy of the PC is effectively reduced to zero 
(Fig.EJ. This elimination of the Coulomb self-energy is essentially due to correlation of PCs 
in the complex and can not be described by Poisson-Boltzmann mean field approximation. 
We define /ic < (c stands for "correlation") as the chemical potential related to the elim- 
ination of the Coulomb self-energy of PC in the complex. Related to the charge of PC, it 
acts as an external voltage overcharging PA. With increasing x, the inverted charge of the 
complex increases. At certain critical x = Xm {itl stands for "maximum charge inversion") 
(see Fig. 121), the maximum charge inversion is achieved where /ic is balanced by the Coulomb 
repulsive energy of the complex to a PC. The topological structure of the phase diagram is 
almost symmetric about x = 1. The only asymmetry appears at a; > Xm- Here additional 
PCs are not attracted to maximum charge-inverted PAs and stay free in the solution (see 
Fig. ED. 

Now let us discuss screening by a monovalent salt. This screening effectively cuts off the 
range of the Coulomb interaction at the distance r^. As decreases, first the long-range 
Coulomb repulsion is reduced, while the short-range correlation induced attraction is not 
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FIG. 4: Possible phases and phase coexistence due to disproportionation. Symbols used here are 
explained in Fig. ^ (a) Coexistence of free complexes and a single drop (inter-complex dispropor- 
tionation): PCs disproportionate themselves among PAs so that a portion of PA-PCs complexes 
are neutral and condensed in a macroscopic drop, while the rest of them are stronger charged and 
free, (b) The tadpole phase (intra-complex disproportionation): PCs disproportionate themselves 
within each PA-PCs complex to form a "tadpole" , with a neutral condensed "head" and a charged 
"tail", (c) Coexistence of tadpoles and free complexes, (d) Coexistence of tadpoles and a single 
drop. In these two cases, inter and intra disproportionation are combined to achieve a lower free 
energy. 

affected. Accordingly, all the ranges of condensation in the phase diagram become wider 
(Fig. 121). Eventually, in the limit of very small r^, the short-range correlational attraction is 
also screened out and the macroscopic drop completely dissolves (not shown in Fig. ^ . In the 
intermediate range of we are interested in, there are two major effects of screening. First, 
the tadpole configuration disappears at certain (Fig. |21). Indeed, if we compare intra and 
inter disproportionation, the former has a lower Coulomb energy but a higher surface energy. 
At large r^, the Coulomb energy is more important and tadpoles are preferred, and vice vera. 
In the analogy to the phase diagram of water, it is like the solid-vapor phase transition at 
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FIG. 5: An illustration of charge inversion of a PA molecule by flexible PCs when they are both 
strongly charged. Negative PA charges are shown by black dots. Positive PC charges are shown 
by white dots. When a new PC molecule is adsorbed to a neutral PA-PCs complex, its charge is 
fractionalized in mono-charges and its Coulomb self-energy is eliminated by redistribution of all 
PCs in the complex. In reality, the numbers of charges of PA and PC can be much larger than 
numbers shown here. Imagine for example that PA and PC have charges — lOOOe and -|-100e. 

low temperature. Second, the single drop phase occupies a finite range of x around x = 1 
at a finite r^, which grows with (Fig- I2))- Recall that at —>■ oo, the macroscopic 
drop should be neutral and therefore it exists only at x = 1. If each condensed complex 
were charged {x ^ 1), the total charge of the macroscopic drop would be proportional to 
its volume and the Coulomb energy per unit volume would be huge proportional to its 
surface. On the other hand, at much smaller than the size of the macroscopic drop, the 
Coulomb energy is not accumulative but additive for each volume of the size r^. Therefore 
the macroscopic drop can tolerate some charge density and the range of the single drop 
phase in the phase diagram widens. 

Currently our theory can be compared with experiments only qualitatively since in many 
cases it is not clear, whether the equilibrium state of the system is reached in experimental 
times due to the slow kinetics. Also the interesting tadpole phase could be very hard to 
realize due to a very large critical Vg above which tadpoles can exist. In solutions of DNA 
with PC, charge inversion of complexes is observed at x > 1 . The size of condensed 
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particles reaches maximum close to x = 1 corresponding to the single drop phase in our 
phase diagram. When x 7^ 1, the size of condensed particles decreases in agreement with 
our equilibrium phase diagram In solutions of DNA with basic polypeptides, at 

X < 1, it is observed that DNA molecules exist simultaneously in two distinct conformations, 
i.e., elongated conformation and condensed conformation [4]. This corresponds to a phase 
coexistence of free complexes and a single drop in our phase diagram (xi < x < 1 in Fig. |2I). 
On the other hand, the enhancement of condensation with the help of simple salt is observed 
in Ref. Certain tadpole-like phases have also been observed in experiments 
although we do not think they are equilibrium tadpole phases discussed in this paper. 

This paper is organized as follows. In SecHIl we discuss all possible phases when ^ 00. 
We then consider the role of screening by monovalent salt and get the complete phase 
diagram in Sec. IIIIl In Sec. IIVI we discuss x > 1 side more carefully and reveal another 
possibility for the phase diagram. In Sec. El we estimate parameters e and /Xc microscopically 
in the case of strongly charged PA and PC. In Se c. IVIL w e discuss the role of translational 
entropy of PC in connection with previous works [l^llJ]. We conclude in Sec. IVIII 



II. PHASES IN THE ABSENCE OF MONOVALENT SALT 

In this section we discuss all possible phases in the absence of monovalent salt, i.e., — >■ 
00. Although this situation is not realistic, it serves as a start point for more complicated 
theory at finite Vg (see Sec. IIII|) . Here and below. We focus on the role of the Coulomb 
interaction and neglect all other interactions such as hydrophobic force. We focus on x < 1 
side and postpone the discussion of x > 1 side till Sec. IIVI 

Let us first clarify the terminology "phase" . By definition, a pha se of a system containing 
one or several substances is an equilibrium homogeneous state [l^. The chemical potential 
of each substance in each phase is uniquely defined and calculated independently. In the 
present case, we have a system of two substances: PA and PC. Following the definition, there 
are three possible phases: the phase of a single drop (Fig. [TJi) , the phase of free complexes 
(Fig.HjD), and the tadpole phase (Fig. |2t). Although we cannot proof rigidly that only these 
three phases can appear in the phase diagram, we have a quite convincing argument. We 
show in the last subsection of this section that three possible phases or phase coexistence 
are ruled out by free energy considerations. 
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At X = 1, there is no Coulomb repulsion between complexes, the phase of a single drop is 
preferred due to the short range correlation energy gain (see Fig. E)). At a; ^ 1, the phase of 
free complexes is preferred because of strong Coulomb repulsion between complexes. Then 
it is natural to conclude that the tadpole phase, which is kind of combination of the other 
two phases, appears in vicinity of x = 1. To testify this idea, we focus on possible phase 
coexistent regions. For example, the phase of free complexes can coexist with the tadpole 
phase, and the tadpole phase can coexist with the phase of a single drop. By minimizing 
the free energies in these two regions respectively, we can find the ranges of x in which one 
or the other phase is preferred, or the two phases coexist. This leads to the phase diagram 
at — >• oo. 



A. The free energies of the three phases 

For pedagogical reason, let us start from the free energy of the tree phases. We first 
notice that the Coulomb energy of each charged complex can be calculated by considering 
it as a conductor with certain capacitance. To see this, let us recall that the total chemical 
potential of a PC molecule adsorbed in PA is given by 

/i = /Xc + #• (1) 

Here the first term is the correlational chemical potential, and the second term is the electric 
energy of a PC given by the local electric potential in the complex. Since both /x and /ic 
are the same along the complex, must be the same in the complex. It is in this sense that 
the PA-PCs complex can be considered as a conductor and the concept of capacitance can 
be used to calculate its Coulomb energy. 

We denote N as the total number of PA molecules in the solution. For the phase of free 
complexes (see Fig. [TJd), we have 



N 



2C 



+ nEin) 



(2) 



Here C is the capacitance of a free complex, —Q and q are the bare charges of PA and PC, 
n is the number of PC in each free PA-PCs complex, and E{n) < is the correlation energy 
of a PC in a complex as a function of n. In this expression, the first term is the Coulomb 
self-energy of free complexes. The second term is the negative correlation energy of PCs 
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in free complexes. We emphasize again that we study very long and strongly charged PC 
and PA such that their translational entropies are negligible. Since all PCs are adsorbed to 
PAs, the net charge of each free complex, {nq — Q), is equal to (x — 1)Q. We will show in 
Sec. |V]that near the phase coexistence region, \n — ni\ '^Ui, where rii = Q / q is the number 
of PC in a neutral complex. Consequently, /Uc(n) = d[nE{n)]/dn is approximately equal to 
its value fic n = Ui. Furthermore, it is convenient to consider the average free energy of 
each complex, / = F/N, instead of F. Finally, we rewrite the free energy as 

/ = ^ + i^c + niE{ni), (3) 

where C is given by 

as the capacitance of a cylindrical capacitor. Here L is the length of a free PA-PC complex 
and -D = 80 is the dielectric constant of water. 

Using similar notations, for a phase of a single drop (see Fig. Ei), we have 

f = n,E{n,)+e. (5) 

Here remember that each complex in the drop is neutral. The expression includes the 
negative correlation energy of PCs in neutral complexes, and the negative correlation energy 
of a neutral complex due to condensation, e. 

For a tadpole phase (see Fig. ^), we denote z as the fraction of the tail part of each 
complex. We have 

, (x-lW , / 36\ (x-l)Q ^, , 

/(^) = ^ + (1 - ( ^ " ) + g /^^ + ^^^(^')- (6) 

This free energy is almost a combination of the last two free energies except two differences. 
First, the capacitance of a tadpole is determined by its tail, given by 

C - (7) 
^* " 21n(^L/6)- 

Second, the correlation energy due to condensation is gained only in the "head" of the 
tadpole. Therefore there is an additional factor 1 — z in the second term. The surface 
energy of the head is also included in the second term R is the radius of the head of 
the tadpole (see Fig. HId), given by 



R 



1/3 



(8) 
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B. The phase diagram at 



oo 



Now we are ready to discuss the free energy in the phase coexistence regions and get 
the phase diagram at oo. We first consider the coexistence of the phase of free 

complexes and the tadpole phase (see Fig. Ht)- We denote y as the fraction of tadpoles in 
the coexistence. Combining Eqs. Q and we have 

(a;-l)2g2 X / 36\ {x - 1)Q 



+ y{l-z)e 1 



-He + riiEirii 



(9) 



2[yCt + {l-y)C] ■ \' Rj q 

Here the first term is the Coulomb energy of the system of free complexes and tadpoles 
with capacitances C and Ct correspondingly. These capacitances are additive because all 
complexes are in equilibrium with respect of exchange of PC and the electric potential of all 



complexes is the same (see discussion at the beginning of the subsection III A|l . 

To Minimize this free energy with respect of two variables y and z, it is convenient to 
separate the free energies into two parts 

1 — x)'^yz In z 5y{l — z^^^ 



where 



1-y + yz 



a 



+ 



a- 



1 -y + yzf 



+ 



DL\e\ 



(10) 



(11) 



is in order of 1 for strongly charged PA and PC (see Sec. |V|), and Eqs. (jD), ((Tj) and (jH)) have 
been used. In this expression, the last two terms of Eq. © have been neglected since they 
are independent of y, z. The first square bracket contains the main terms of the Coulomb 
energy and the correlation energy, in which In^L in Ct has been replaced by InL and the 
surface energy has been ignored. The second square bracket contains correction terms in 
the first order to make up the two neglects in the main terms. We will see below that the 
small parameter we used to expand this free energy is \ii{L/b)/{L/bY^^ , for L ^ b. 
We first minimize the main terms in Eq. (|TO|l and get 



y{l - z) = 1 - {1 - x)^a\n{L/b). 
Putting it back to the correction terms in Eq. ()10|) and taking the minimum, we get 



(12) 



z = 1 



"51n(L/6)" 
_3(L/6)i/3_ 



3/4 



(13) 
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Combining the above two conditions, we have 



y 



1 - {1- x)Ja\n{L/b) 



3{L/b) 



1/3 



3/4 



(14) 



51n(L/6) 

Specifically, considering y = and y = 1, we get the boundaries of the phase coexistent 
region (see Fig. ^ 

1 



1 — Xi{00) 
1 - X2(oo) 



aln{L/b) 
1 



51n(L/&) l ^/'^ 



(15) 
(16) 



'a\n{L/b) 

Here and below x(oo) means x{rs oo). When x increases from xi(oo) to ^2(00), y 
increases from to 1 linearly. This is actually the level rule and the sign of a first order 
phase transition. 

For a very long PA, L ^ b, parameter \n{L/b)/ {L/bY^'^ <^ 1. Therefore z is close to 1, 
the free energy expansion in Eq lfTIH) is self-consistent. Also a;i(oo) is close to 1 and X2{oo) 
is close to xi(oo). 

We consider the other phase coexistent region (see Fig. |lji) following the same approach. 
Combining Eqs. ^ and (jHl), we have 

f{y, z) = ^^^^ + yil - z)e f 1 - |U (1 - y)e + + n,E{n,). (17) 



The free energy similar to Eq. (fTUI) is 

[l-xf\xi{L/b) 



f{y.z)/\e\ 



a- 



yz 



+ yz-l 



+ 



'^{l-x)^\nz ^ 5?/(l- 2)2/3 



yz 



{L/byr^ 



(18) 



Minimizing the main terms (in the first square bracket), we get 



yz = [1 — x)Ja\n{L/b). 



(19) 



Putting it back to the correction terms (in the second square bracket) and taking the mini- 
mum, we get 



51n(L/6) 
(L/6)i/3 • 



Combining the two conditions, we have 



:i-x)y^(L/6)i/3 



5Jln{L/b) 



(20) 



(21) 
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Therefore the boundaries of this coexistent region are (see Fig. 12)) 

, , 1 51n(L/6) , , 

1-0:4(00) = 0. (23) 

We see that the same small parameter ln(L/6)/(L/6)^/'^ works here. Again, the phase 
transition from the tadpole phase to the phase of a single drop is in the first order. 



C. Why not other phases 



In this subsection, we argue that several other possible phases or phase coexistence do 
not appear in the phase diagram at ^ 00. 

One possibility would be that the tadpole phase does not appear at all, but the coexis- 
tence of the other two phases shows up in vicinity of x = 1. This is called inter-complex 
disproportionation (Fig. |3^). We showed in Ref. 15| that this coexistence has a higher free 
energy than the tadpole phase in some range of x < 1. As far as the tadpole phase exist, 
it immediately follows that the only phase coexistent regions can exist are those involving 
tadpoles (see Fig. Et,d). And the inter-complex disproportionation region (Fig. 0^) has 
nowhere to show up. 



(a) 



(b) b 




FIG. 6: Comparison of two- heads and one-head configurations, (a): a two-head configuration. The 
diameter of the smaller head is d. (b): a one-head configuration made from (a) by combining the 
two heads then releasing tail of length order of d from the head, such that the two configurations 
have same capacitance. The Coulomb energy is the same for the two configurations, but the surface 
energy is higher in (a). 

The second possibility is a phase in which every complex has more than one "head" 
and therefore forms a necklace structure. Let us argue that the one-head-one-tail tadpole 
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configuration is the best for intra-complex disproportionation. To do tliis, we liave to include 
tfie capacitances of lieads which gives a small correction to Eq. ((7j) . First of all, let us show 
that for a complex with given charge, one head is better than two heads. Consider an 
arbitrary two-heads configuration with the diameter of the smaller head d {Fig^). We 
can always construct an one-head configuration from it by combining the two heads then 
releasing additional tail of the order of d from the head in such a way that the capacitances 
of the two configurations are the same (Fig|ni3). The total free energy consists of the long- 
range Coulomb energy and the short-range correlation energy in droplets. For the two 
configurations, the Coulomb energies are the same since the capacitances are equal. But 
the surface energy is higher in the two-heads configuration since the surface is much larger. 
Thus, for any two-heads configuration, we can always find a one-head configuration with 
lower energy. By similar argument, obviously a configuration with many heads along the 
complex is even worse. Furthermore, the single head always prefers to be at the end of the 
tail. This can be understood by considering a metallic stick with fixed charge on it. The 
electric field is larger at the end of the stick than in the middle. Therefore to reduce the 
Coulomb energy, it is better to put a metallic sphere at the end of the stick to make field 
there smaller. 




FIG. 7: An illustration of instability of hypothetic micelle-like droplets in a PA-PC system. Symbols 
are explained in Fig. ^ If tadpoles are not stable, one could think that they merge into micelles. 
But micelles would immediately break down to reduce Coulomb energy. 

The third possibility is a phase in which every few complexes come together to form a 
micelle-like object (see Fig. [7j). However, micelle-like droplets are not stable. For example, 
suppose several charged PA-PCs tadpoles would gain energy merging into a micelle. Instead 
of sharing one head by several tails, to reduce the Coulomb energy, it is better to break 
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the micelle into one tadpole and several tails. Repeating this process, certainly one ends up 
with a two phase coexistence (see Fig. [7j). 

Finally, we want to point out that the necklace phase and the micelle phase are always 
worse in free energies than the tadpole phase. This is not only true at — > oo. One can 
show that they do not exist in the case of finite too. 



III. PHASE DIAGRAM IN A PLANE OF SCREENING RADIUS AND CHARGE 
RATIO 

In this section, we consider a more realistic situation of finite and get the complete 
phase diagram in a plane of h/rs and x (see Fig. Again we focus on x < 1 side. 
The Tg range we are interested in is 6 <^ ^ L. As becomes smaller than L, the 
long-range Coulomb repulsion is first reduced while the short-range correlational attraction 
is not affected. Eventually, when < 6, the short-range correlational attraction is also 
screened out and the macroscopic drop completely dissolves but we do not consider such salt 
concentrations. In the range of our interest, 6 ^ ^ L, there are two major implications 
of screening (see Fig. 12)). First, the phase of single drop grows up. Second, the tadpole phase 
and related phase coexistent regions are destroyed. 

Let us first discuss the general confinement to the phase diagram given by the Gibbs' 
phase rule Suppose the number of coexistent phases is m. Then for our two substances 
system, the number of independent equations following the condition of equal chemical 
potentials is 2(m — 1). The number of unknowns is m + 1. For example, they can be chosen 
as the charge ratio x in each phase and the common parameter fe/r^ shared by all coexistent 
phases. To have a solution to the equations, we require 2(m— 1) < m + 1. This gives m < 3. 
That is to say, the number of coexistent phases cannot be more than 3. It is convenient to 
define a chemical potential 

A. = 5^. (24) 

and consider the phase diagram on a plane of and fo/r^ (see Fig. |H)). When three phases 
coexist, all variables are completely determined. It corresponds to a point in the phase 
diagram (the triple point). When two phases coexist, there is one thermodynamic degree of 
freedom. It corresponds to a line in the phase diagram. If we change variable from /i^ to its 
conjugate variable, x, we get our phase diagram (Fig. 12)). One can easily show that in this 
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phase diagram, three phases coexist on a hne, and two phases coexist in a region. 



1^ 



X 




b 



FIG. 8: Schematic phase diagram on a plane of fix and r^. Here fix is a chemical potential which 
is conjugate to x. 

The phase diagrams we get is very much similar to the phase diagram of water Q|. 
Indeed, if we consider fix as pressure, r<./6 as temperature, and x as volume, then Fig. |Slis 
like the FT diagram, and Fig. |21is like the TV diagram of water. Notice that in our system, 
real temperature or pressure plays no role. The system is always at room temperature. 
And the pressure is negligibly small, since in comparison of the Coulomb energy, we have 
completely ignored the translational entropies of PA and PC, which is related to the volume 
and the pressure of the system. 

Having figured out the topological structure of the phase diagram, let us determine all 
phase boundaries by considering the three two-phase coexistent regions. We start from the 
coexistence of tadpoles and free complexes (Fig. |^). In the presence of monovalent salt, the 
free energy is still given by Eq. Q, but the expressions of C and Ct can be different. Let 
us first assume that zL < Tg <^ L. Then Ct is still given by Eq. ((7j), while C is now 



21n(r,/6)' 

Following the same procedure as in subsection 111 Bl we get 

Ts 51n(rs/6) 




(25) 



2.7l''''P[ 3(L/6)i/3 



(26) 




(27) 
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At y = and y = 1, we get the boundaries of this coexistent region (see Fig. ^ 

l-x,(r,) = -=i=, (28) 



l-a;2(r-s) = , " exp 

2.7LJa ln{rjb) 



5 ln(rs/6) 
'3(L/6)V3 



(29) 



In the hmiting case of — L, these equations go back to Eqs. (fT^ and (fT^ . 

If Ts <C zL, not only C, but also Cj gets a new expression 

^* = 2ln{rjby 

In this case, one can check that the free energy given by Eq. does not have a minimal 
extremum. Indeed, the Coulomb energy of the system becomes so short-ranged that there 
is almost no difference to distribute all charges to free complexes or tadpoles. The surface 
energy of heads dominates and the coexistence of free complexes and the single drop is 
preferred (see Fig. |2I). 

We now consider the coexistence of tadpoles with the single drop (Fig. HJi). In the case 
of finite r^, the free energy given by Eq. (|T7jl should be revised to 

jx-lYQ' ( 3b\ , (x-l)g 

f{y, z) = ^^yC, + [l-y)C'] + -'^'y-R)^ -^^^ + 

Here C represents the capacitance of a condensed complex in the single drop. As discussed 
in SecUl when ^ oo, the macroscopic condensate is almost neutral and the phase of single 
drop exist only at x = 1. For finite r^, the macroscopic drop can tolerate certain charge 
density. Therefore C appears. In order to calculate C", we assume that the macroscopic 
drop is u„,for,„.y charged 0. If the charge dcrs.ty of the macroscopic drop ,s , and the 
charge of each complex is i^b Lp/A the electrical potential of the macroscopic drop is 

°° P^Anr^dr = (32) 
Q Dr D ' ^ ' 

This gives 

_ T^b^Lp _ Db'L 
^ " 40 " 16r2 • ^'^'^^ 

When ^ oo, the capacitance C" ^ as expected. 

In the case of zL < Ts L, Ct is still given by Eq. (jZj). Minimizing the free energy as in 



subsection IIIBl we get 

^ 51n(L/&) ln^(L/&) 

(L/6)i/3 + 8(r,/6)2 • ^ 
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And y is still given by Eq. (fT^. Setting ?/ = and y = 1, the boundaries of the coexistent 
region are (see Fig. 12)) 



51n(L/6) 



1 - x^ivs) = 0. 



1 + 



(L/6)i/=^ ln(L/6) 
40(r,/6)2 



(35) 
(36) 



In the other case of ^ zL, Ct is given by Eq. ()30p. Again the free energy Eq. ()31|) has 
no minimal extremum but a maximal extremum. And the minimal free energy is achieved 
at 2 = 1. This leads to the coexistence of free complexes and the single drop is preferred 
(see Fig. EI). 

We are now ready to determine the position of the line at which three phases coexist in 
Fig. 121 (it is also the upper boundary of the tadpole phase and the two two-phase coexistent 
regions). The critical r^, Tc, is determined by the cross point of X2{rs) and x^lrs). In the 
leading order, we have 

r.c:^ L^/%'^Hn{L/b). (37) 

Unfortunately, Vc is large and not easily realized in experiments. For example, for b = 7 A and 
L = 1006, Tc — 996 = 693 A. This means that the tadpole phase and all phase coexistence 
related with it are typically not relevant to experiments. Physically, comparing with inter- 
complex disproportionation, the tadpole phase (intral-complex disproportionation) is better 
in the Coulomb energy, and worse in the surface energy. For the Coulomb energy to be more 
important, has to be large. Therefore the tadpoles must be extinct at small Vg. 

Finally we consider the coexistence of free complexes and the single drop (see Fig. |3K) 
when rg < Vc- Combining Eqs. and (jS)), we have 

(x 



f{y) 



' - + {l-y)e+ ^ -^/i, + riiEirii 



2[yC+{l-y)C'\ q 
where C is given by Eq. ()25|). Minimizing the free energy, we have 



{x - If 



ln(rs/6) 8rl 



D\6\L 



ln(r^/6) 



+ 



;i - y)b^ 



And 



1 - xJTs < Tc 



ln{rs/b) 



/ln(r,/6) 



a 



(38) 

(39) 

(40) 
(41) 
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Here we put < Tc as the argument to remind us that these expressions are meaningful 
only at < re- We see that the width of the single drop phase grows proportionally to 
with decreasing (see Fig. |2)). When — > rc, X4{rs) ^ 1 as expected. 



IV. POSSIBLE NEW PHASES AT x > 1 SIDE 



In this section, we discuss the phase diagram at x > 1 side. At a; > 1, the number of 
PC molecules is more than enough to neutralize all PA molecules. The signs of charges of 
free complexes or the tails of tadpoles are inverted to positive. In spite of this difference, all 
physics we discussed in the last two sections remain valid. Therefore one expects that the 
topological structure of the phase diagram is symmetric about x = 1. And the boundaries 
x[{rs) (see Fig. 12]) satisfy 

x[-l = l-x„ (42) 

where i = 1,2, 3, 4. 

However, a little asymmetry exists since the charge inversion process cannot go forever, 
but reaches its maximum value at certain critical value of x, Xm- Xm is determined from the 
balance of the gain in the correlation energy with the overall Coulomb repulsive energy for 
a PC molecule. For example, for a free complex, we have 

I^J = '^"-"Q" . (43) 

where each free complex carries net charge {xm — 1)Q- In this equation, the left hand side 
is the magnitude of the gain in correlation energy, and the right hand side is the Coulomb 
repulsive energy given by the net charge of the complex. At a; > Xm, extra PC molecules 
are not attracted to the complexes, but stay free in solution (see Fig. |21). Using Eq. (j^ . we 
get 

When Ts ^ CO, 

'''' 

Furthermore, the existence of Xm gives another possibility to the phase diagram. For 
our purpose, it is enough to consider the simple case where the free complexes and the 
single drop coexist (see Fig. One can show that the conclusion is the same for phase 
coexistence involving tadpoles. 
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Indeed, what we discussed above is self-consistent if x'^ < Xm- But what if x[ > Xm? To 
answer this question, we first notice that x[ has another physical meaning. In Eq. (jH^ . for 
the accuracy we needed, keeping the first term on each side, we get the net charge of each 
free complex in the phase coexistent region 

(^-1)^ = Q (46) 

y ^aln{rs/b) 

According to Eq. (pOj). this net charge is equal to {x[ — l)Q. When x[ > Xm, to keep the 
phase coexistence, each free complex should carry charge {x'l — 1)Q while it is only allowed 
to carry {xm — 1)<5 because of the finite correlation chemical potential /ic- In this situation, 
x[ loses its physical meaning and the charge of each free complex saturates at the maximum 
value {xm — The free energy of Eq. (jnH|) should be revised to 

m = 2/^^^^^^^ + + (1 - y)e + nMn.). (47) 

Here for the purpose of discussion, the second order term related to C has been ignored. 
With help of Eqs. pn|) and this free energy can be written as 

f{y) = -y ^ ^ - (1 - y y ' J ^ + Nn,E{n,). (48) 

Clearly, when x[ > x^, the minimum of f{y) is reached at y = 0. Therefore at x > 1 side, 
we arrive at a phase of total condensation in which all complexes are condensed but some 
PCs are free. This leads to a different phase diagram (Fig. E]). 

It is convenient to classify all PA-PC systems into two categories. In the first category, 
x'^(oo) < Xm(oo), and we have phase diagram Fig. El In the second category, x[{oo) > 
Xm(oo), and we have phase diagram Fig. El Interestingly, both two categories are realistic 
as we will see in Sec. El Notice that in the second case, decreasing to a critical value. 
To, eventually leads to inversion of inequality x[ > Xm (see Fig. (HI). We discuss this effect in 
detail in the next section after we estimate e and fic microscopically. 

In the case when x'i{oo) > Xm{oo), the phase of single drop expands around x = 1 with 
growing L/vg ai x > 1 side too (see Fig. [HI). We calculate the boundary of this phase with 
the phase of total condensation at > tq. What we need to find out is just how many 
excessive PCs the macroscopic drop can tolerate at finite r^. Applying the condition of 
maximum charge inversion to a condensed complex, similarly to Eq. (021), we have 

|,,| = (49) 
21 




FIG. 9: Phase diagram of solution of PA and PC in the case when at — > oo, we formally have 
x'l > Xm- The meaning of axes and symbols are the same as Fig. |^ 



x[{r,<ro) = l + ^Bhr- (50) 



Therefore 

D\iic\h'^L 

At Ts < To, we arrive at inequality x'i{r-s) < Xm{rs) and Eq. (jlUj) gives the boundary of the 
phase of single drop (see Fig. We will discuss the transition at = ro in detail in the 
next section. 



V. PHASE DIAGRAM OF STRONGLY CHARGED POLYELECTROLYTES 

In this section, we consider a simple system where linear charge densities of PA's and 
PC's are equal. Both of them are strongly charged such that every monomer carries a 
fundamental charge e and e^/Db ~ ksT (see Fig. EJ. We estimate parameters e and /ic 
microscopically and choose from the two phase diagrams shown in Figs. |21 and IHl 

We first consider the case of Ts — > oo. As discussed in Sec. HI fic is equal to the Coulomb 
self-energy of a PC, 

/^c = -|^ln(g/e). (51) 
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On the other hand, when neutral PA-PCs complexes condense, they form a strongly corre- 
lated liquid. Monomers of two PEs locally form NaCl-like structure such that the energy in 

I — I 

order of —e^/Dh per monomer is gained (Fig. [18]. Therefore 



Qe 

e ^ . 

Dh 

We see that indeed a defined by Eq. (fTT|) is in order of 1. 
Substituting L = Qb/e into Eqs. (fT^ and (ji^ . we have 



(52) 



a;;(oo) ~ 1+ , ^ =, (53) 



Accordingly, we get the critical value of q at which 

ge = exp (2^HQje)^ . (55) 

When q > qc, or PC is long enough, x'l < Xm, we have the phase diagram of reentrant 
condensation (horizontal axis of Fig. El). When q < q^, or PC is relatively short, x'l > Xm, 
we have the asymmetric phase diagram with total condensation at a; > 1 side (horizontal axis 
of Fig. in)). The possibility of having two different phase diagrams for different q is related 
to the interplay of the logarithms in Eqs. (j5T| and (j5l|l . For fixed Q, when q increases, Xm 
increases but x[ is fixed. Therefore we can have either Xm < x'l or Xm > x'l by changing q. 
Notice that qc is exponentially smaller than Q. 

Now let us consider the effect of screening by monovalent salt. We are interested in the 
case of Tg ^ b when the short-range correlation is not affected yet and e is fixed. According 
to Eq. 

x[{rs)^l+ , ^ (56) 



In order to discuss Xm(rs), we consider two different cases, qb/e <^ <^ Qb/e and 
& ^ ^ qb/e. When qb/e <^ ^ Qb/e, the chemical potential fic is still given by 
Eq. (jni)- From Eq. 

Accordingly, the critical value of q at which x[{rs) = Xm{rs) is 

qdrs) = exp \2J\n{r,/b)\ . (58) 
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It decreases with decreasing r^. Therefore, for a system with q < gc(oo), q can be larger than 
qd^s) at small Vs- Correspondingly, at x > 1, for small r^, the phase of total condensation 
is replaced by the phase coexistence of free complexes and the single drop. We have a phase 
diagram shown in Fig. 13 Letting q = qdrg), we get the critical value of at which this 
phase transition happens 



Notice that ro is much larger than qb/e. 

In Fig. El x[{rs) > Xm{rs) at > Tq, while x[{rs) < Xmi^s) at r., < {x[{rs) curve at 
Ts > ro is not shown). By definition of Tq, the curves x[{rs) and Xmi^s) merge at = rg. 
The curve xd^s) is given by Eq. (j5(J|) for > Tq and Eq. (j4(Jj) for < Vq. At = Tq, 
these two expressions are equal to each other. At a; > 1 side, the solid line at L/vg = L/tq 
corresponds to a first order phase transition. Notice that tq can be either smaller or larger 
than Tc- Here for simplicity, only the former case is shown in Fig. El 

When h <^rs <^ Qb/e, 



and Xm{rs) = 3/2 |l2|. In this case, we always have x'^Tg) < Xm{rs) as shown in Figs. El 
andU 

Finally, in all cases discussed above, the value of x'^ is close to 1, i.e., \n — ni\ Ui in 
the condensation regime. Therefore the approximation used in Sec. IHlthat fic is a constant 
is valid. 

VI. THE ROLE OF THE TRANSLATIONAL ENTROPY OF POLYCATIONS 

A major approximation in this paper is that the translational entropy of PCs is negligible 
(we can always ignore PA's translational entropy since it is much longer than PC). In this 
section, we would like to discuss the validity of this approximation and the role of the 
translational entropy. 

First, let us estimate when this approximation is valid. Consider PCs with concentration 
p in the solution. The free energy due to its translational entropy is ksTlnlpvo), where Vq is 
the normalizing volume. One the other hand, according to Eq. (|5ip. the Coulomb energy is 
in the order of —qe/Db ~ —qkBT/e. They are equal at the critical value, p = exp(— g/e)/fo. 




(59) 



qe 
Db 



ln(rs/6) 



(60) 
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Therefore for a long PC with large g, we can ignore PC's translational entropy even at 
exponentially small p. 

If PC is very short, its translational entropy should be included. DNA with short 
polyamines is a good case, the phase dia- 

gram gets another dimension, say, the concentration of PC, p. The effect of PC entropy was 
discussed in detail in Ref. (l3l . Il4| in which the phase diagram is drawn in a plane of two 
concentrations of oppositely charged colloids at given r<j. Here we discuss the same effect 
in the language of total charge ratio x used in this paper in the simple case where x[ > Xm 
and Ts — oo. For simplicity, we neglect the possibility of intra-molecule disproportionation 
and the tadpole phase. 

In this case, the free energy in Eq. (j38|) gets an additional term due to the translational 
entropy of PCs Q| 

'{nq-Qf 



fin,y) = y 
+ 



2C 



+ nEin) 



xQ 



+ {I - y)[niE{n{) + e] 



In 



' ynq 
P 



X 



-p 



(61) 



where e is the natural exponential. Here the expression in the square bracket before the 
logarithm represents the number of free PC in the solution, while the expression in the round 
bracket in the logarithm represents their concentration. 

Now n and y are two independent variables. Taking dF/dn = and dF/dy = 0, we get 



In 



ynq 
p —p 
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xQ' 
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-P Vq 



{nq - Qf 
2C 



+ {n- rii) <^ /Xc - In 



{nq - Q)q 
C ' 
/ ynq 



1-!/ 



X 



-P vo 



(62) 
(63) 



In Eq. (j62|) . eliminating n by Eq. (j63|) . we can calculate the boundaries of the condensation 
regime by setting y = and y = I. For ?/ = 0, we get two boundaries of a total condensation 
phase, 



p[l-- 

X 

1 

X 



pll 



Pi, 
p'l- 



For ?/ = 1, we get two boundaries of the free complexes phase, 

p{^-"i) =Pu 



(64) 
(65) 

(66) 
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FIG. 10: Phase diagram in a plane of PC's concentration p and total charge ratio x {x'l > Xm and 
rg oo). Symbols used are explained in Fig. ^ It shows how the phase of total condensation 
replaces that of partial condensation with decreasing p. 



Here 
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(67) 

(68) 
(69) 



Accordingly, as shown in the phase diagram Fig. ^[ a regime of total condensation is 
sandwiched by two regimes of partial condensation, which are further sandwiched by two 
regimes of free complexes. 

In Fig. ^1 results of previous sections are recovered in the limit p — >■ oo. Due to the 
translational entropy of PC, at finite p all critical x are shifted to higher values. At the same 
time, a total condensation region acquires a finite width even in the absence of monovalent 
salt 13]. In the limiting case where x — > oo, the concentration of PA is much smaller than 
that of PC, and the entropy of PC is fixed, which offers a fixed charging voltage to PA. 
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As a result, all PA-PCs complexes are either totally condensed or totally free. The partial 
condensation regime disappears jisl - 

VII. CONCLUSION 

In this paper, we discussed complexation and condensation of PA with PC in a salty water 
solution. Using ideas of disproportionation of PCs among complexes and inside complexes 
(inter- and intra-complex disproportionations) we arrived at the two phase diagrams in a 
plane of x (ratio of total charges of PC and PA) and L/rg (ratio of the length of PA, L, 
and the Debye-Hiickel screening radius, Vg) shown in Figs. ElandlHl In the case of strongly 
charged PA and PC, we find that both two phase diagrams are possible depending on the 
relative length of PC to PA. Fig. [21 corresponds to a more generic case of relatively long 
PC, while Fig. El to the case of relatively short one. Our phase diagrams show how total 
condensation is replaced by the partial one and then by phases of stable complexes when x 
moves away from x = 1. 

We discovered two new features of the phase diagrams. First, at large screening radius 
they include a new phase of tadpoles and corresponding phase coexistence. Second, we 
found that the phase of the single drop formed at x close to 1 widens with decreasing as 

Although we talked about strongly charged PA and PC one can also consider phase 
diagram of weakly charged PA and PC and develop a microscopic theory for it. In both 
cases, the qualitative picture is the same since our discussion of the phase diagram is rather 
general and independent of the microscopic mechanism of the short-range attraction. 

As mentioned in the introduction, the problem we solved in this paper should be con- 
sidered as an example of a more general problem of the phase diagram of the solution of 
two oppositely charged colloids. Another important system of this kind is a long PA with 
many strongly charged positive spheres. When long double helix DNA plays the role of PA, 
this system is a model for the natural chromatin. Therefore, we call such system artificial 
chromatin Our phase diagrams with all new features including tadpoles should be valid 
for artificial chromatin as well. 

There is another class of systems where only some of our predictions are applicable. In 



the Ref. 



1J| we considered solution of large negative spheres with positive spheres which 
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are smaller in both radius and charge. Complexation and condensation of such spheres obey 
the phase diagrams similar to discussed above. For example, screening by monovalent salt 
again leads to expansion of the range of the single drop (total condensation) phase. 
Another our prediction, the tadpole phase, however, is not applicable to this case, because 
it is essentially based on the polymer nature of PA. 

In a case when the role of PA is played by DNA, one should remember that the double 
helix DNA is so strongly charged that the effect of the Manning condensation by monovalent 



counterions must be included |13|. Since DNA complexes with positive macroions (PC, 
positive spheres or multivalent cations) , this Manning condensation can be weaker than free 
DNA due to counterion releasing. The quantitative description of this effect depends on 
the geometry of the positive macroions and the microscopic structure of the complex in the 
system jisl]. Generally speaking, this effect leads to the renormalization of the bare charge 
of PA, Q. In this case, total condensation still happens around x = 1 but renormalized 
charge enters in calculation of x. This means that if on the other hand x is evaluated using 
the bare charge of DNA, all phase diagrams are centered around a smaller than 1 value of 

X. 

Our phase diagrams deal with equilibrium states of the system. But not all of them 
can be achieved in experimental time scale due to slow kinetics. Therefore it is not easy 
to directly compare our theory to experiments. For instance, a phase of many condensed 
particles with finite size is often found in experiments which does not appear in our phase 
diagram Q fl fl BLWe believe tl,at tl,. pl,ase is not a real equUibdum state, but the state 
frozen kinetically |2fll . 12 1| . Thus, kinetics is extremely important for applications and we 
plan to address it in the future. 
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